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ABSTRACT 


A  relationship  between  a  person's  lean  body  mass  and  the  amount  of  maximum 
torque  that  can  be  produced  with  each  isolated  joint  of  the  upper  extremity  was 
investigated.  The  maximum  dynamic  isolated  joint  torque  (upper  extremity)  on  14 
subjects  was  collected  using  a  dynamometer  multi-joint  testing  unit.  These  data  were 
reduced  to  a  table  of  coefficients  of  second  degree  polynomials,  computed  using  a 
least  squares  regression  method.  All  the  coefficients  were  then  organized  into  look-up 
tables,  a  compact  and  convenient  storage/retrieval  mechanism  for  the  data  set.  Data 
from  each  joint,  direction  and  velocity,  were  normalized  with  respect  to  that  joint's 
average  and  merged  into  files  (one  for  each  curve  for  a  particular  joint).  Regression 
was  performed  on  each  one  of  these  files  to  derive  a  table  of  normalized  population 
curve  coefficients  for  each  joint  axis  direction  and  velocity.  In  addition,  a  regression 
table  which  included  all  upper  extremity  joints  was  built  which  related  average  torque 
to  lean  body  mass  for  an  individual.  These  two  tables  are  the  basis  of  the  regression 
model  which  allows  the  prediction  of  dynamic  isolated  joint  torques  from  an 
individual’s  lean  body  mass. 


1.0  INTRODUCTION 


A  relationship  exists  between  a  person's  lean  body  mass  and  the  amount  of  maximum 
torque  that  can  be  produced  with  each  isolated  joint  of  the  upper  extremity.  The  use  of 
an  easily  measured  parameter  (lean  body  mass)  to  predict  dynamic  isolated  joint 
torque,  which  is  time  consuming  to  measure,  would  be  extremely  valuable.  In  this 
study,  we  have  collected  maximum  dynamic  isolated  joint  data  on  14  subjects  for  the 
upper  extremity  and  formulated  a  regression  model  which  will  allow  prediction  of  a 
joint  angle  versus  torque  curve  for  a  particular  individual.  There  are  three  phases  in 
our  study:  1)  data  collection,  2)  data  reduction,  and  3)  model  formulation. 

Phase  one  of  our  study,  data  collection,  involved  measuring  the  maximum  torque  for 
all  the  upper  extremity  joints  (shoulder,  elbow,  and  wrist)  at  4  velocities  for  14  subjects. 
Data  for  three  subjects'  lower  extremities  (hip,  knee,  and  ankle)  were  also  collected. 
Torque  was  measured  by  using  a  LIDO  multi-joint  testing  unit  (Loredan  Biomedical, 
Inc.,  West  Sacramento,  California).  The  subjects  were  positioned  so  that  the  axis  of  the 
joint  was  directly  in  line  with  the  axis  of  the  dynamometer  goniometer.  Dynamometer 
attachments  were  selected  and  used  to  isolate  the  joint  being  measured.  In  this 
manner,  all  joints  were  characterized  for  all  axes  of  rotation  over  a  range  of  velocities. 
The  data  were  collected  using  the  Loredan  software,  LIDO  Active  3.3,  on  an  IBM  PC. 
For  all  cases,  the  data  set  consisted  of  torque  and  angle  pairs.  In  addition, 
anthropometric  data  were  also  collected  which  included  height,  weight,  age,  sex,  skin 
fold  measures,  and  dimensional  assessment  according  to  the  format  specified  in 
NASA  Man-Systems  Integration  Standards  (MSIS)  document  [1]. 

The  second  phase  of  the  project,  data  reduction,  began  with  transferring  the  data  to  a 
UNIX-based  workstation  (Silicon  Graphics).  These  data  were  formatted  into  an  ASCII 
file,  noise  filtered,  reformatted,  and  reduced  to  a  table  of  coefficients  of  second  degree 
polynomials.  The  polynomial  coefficients  were  computed  using  a  least  squares 
regression  method.  These  polynomials  represent  the  torque  as  a  function  of  angle 
(i.e.,  torque  =  a+  b*angle  +  c*  angle**2,  where  a,  b,  and  c  are  the  polynomial 
coefficients).  All  the  coefficients  were  then  organized  into  look-up  tables.  These 
tables  represent  a  compact  and  convenient  storage/retrieval  mechanism  for  our  entire 
data  set  and  are  available  upon  request. 

The  third  phase  of  our  project  involved  model  formulation.  Data  from  each  joint, 
direction  and  velocity,  of  an  individual  were  normalized  with  respect  to  that  joint's 
average  and  merged  into  files  (one  for  each  curve  for  a  particular  joint).  Regression 
was  performed  on  each  one  of  these  files  to  derive  a  table  of  normalized  population 
curve  coefficients  for  each  joint  axis,  direction,  and  velocity.  In  addition,  a  regression 
table,  which  included  all  upper  extremity  joints,  was  built  which  related  average  torque 
to  lean  body  mass  for  an  individual.  These  two  tables  are  the  basis  of  the  regression 
model  which  allows  isolated  joint  curve  prediction.  Because  of  the  limited  number  of 
subjects  (3)  for  the  lower  extremity  data  set,  no  correlation  to  lean  body  mass  is 
presented  here  for  these  measurements. 
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Finally,  we  have  encapsulated  the  results  of  this  study  in  a  tool  kit  of  software  routines 
executing  on  a  variety  of  platforms  such  as  UNIX,  DOS,  and  Macintosh  machines 
(Appendix  D).  This  code  contains  our  isolated  joint  torque  model  (with  all  the  tables  of 
torque  coefficients)  and  will  allow  the  prediction  of  dynamic  isolated  joint  torques  from 
an  individual's  lean  body  mass. 
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2.0  OBJECTIVE 


Prediction  equations  have  been  developed  from  lean  body  mass  for  isometric/static 
strength  [2].  Few  studies  have  utilized  lean  body  mass  to  predict  isokinetic  or  dynamic 
strength  [3].  An  absence  of  literature  exists  when  correlating  lean  body  mass  to 
isolated  joint  isokinetic  mean  torque  or  torque  over  an  entire  range  of  joint  motion.  In 
addition,  lean  body  mass  has  not  been  used  to  predict  the  mean  torque  for  dynamic 
complex  tasks  involving  multiple  joints.  The  use  of  an  easily  measured  parameter 
(lean  body  mass)  to  predict  dynamic  isolated  joint  torque  would  be  of  significant 
value. 

The  objective  of  this  project  is  to  develop  prediction  equations  which  can  be  used  to 
calculate  isolated  joint  torque  (either  mean  torque  or  torque  as  a  function  of  angle)  and 
mean  torque  for  a  complex  task  from  the  measurement  of  a  person's  lean  body  mass. 

Specific  aims: 

1 .  Document  all  data  measurement  and  data  processing  techniques. 

2.  Develop  prediction  equations  for: 

a.  Mean  dynamic  torque  over  an  entire  range  of  motion  for  a  particular 
joint  (shoulder,  elbow,  and  wrist)  and  during  a  complex  task  (ratchet 
wrenching  or  an  extravehicular  task). 

b.  Dynamic  torque-position  curves  over  the  entire  joint  range  of  motion 
from  lean  body  mass  for  individual  isolated  joints  (shoulder,  elbow,  and 
wrist). 

3.  Develop  a  set  of  software  tools  which  simplify  dynamic  torque  data  access, 
manipulation,  and  prediction  for  the  set  of  isolated  joints  measured. 
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3.0  METHOD 


3.1  Data  Collection 
Subjects 

Fourteen  subjects  (8  males  and  6  females)  aged  21  to  28  years  volunteered  and 
participated  in  this  investigation.  The  study  was  evaluated  and  approved  by  the 
Institutional  Review  Board  of  the  University  of  Texas  Medical  Branch  at  Galveston, 
Texas.  All  subjects  were  informed  of  potential  risks  and  signed  a  consent  to 
participate  in  the  study. 

Equipment 

The  LIDO  multi-joint  testing  unit  (Loredan  Biomedical,  Inc.,  West  Sacramento, 
California,  Figure  1)  is  an  integrated  system  consisting  of  a  dynamometer  connected  to 
a  personal  computer  via  RS232c  lines.  The  force  unit  comes  with  a  series  of 
attachments  and  a  subject  bench  which  allows  isolation  of  particular  joints.  The 
software  allows  precise  control  of  the  actuator  head  and  the  various  modes  of 
operation  (e.g.,  isometric,  isokinetic,  and  concentric).  In  addition,  a  database  of 
subjects  can  be  maintained  and  displayed  with  the  provided  graphical  software.  The 
unique  feature  of  this  system  is  that  it  outputs  all  data  in  a  machine-readable  form  for 
more  accurate  data  analyses.  This  system  was  used  to  measure  all  the  isolated  joint 
forces  as  well  as  the  composite  test  cases. 

Data  Collection  Procedures 

Data  collection  for  this  project  occurred  over  an  8-week  period  at  the  University  of 
Texas  Medical  Branch,  School  of  Allied  Health  Sciences,  Human  Performance 
Laboratory.  The  general  procedure  for  evaluating  all  the  upper  and  lower  extremity 
joint  movements  and  simulated  ratchet  maneuvers  was  the  same.  Specific  subject 
and  joint  positioning  for  the  isometric  and  isokinetic  tests  is  described  in  each 
independent  section.  The  time  required  to  perform  all  of  the  dynamic  and  isometric 
shoulder,  elbow,  and  wrist  joint  measurements  was  1 .5  to  2  hours;  evaluation  of  the 
hip,  knee,  and  ankle  measures  required  1 .0  to  1 .5  hours;  and  evaluation  of  the 
3  ratchet  wrench  maneuvers  required  0.5  to  1 .0  hours  per  subject. 
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Backrest 


Actuator 


Controller 


Figure  1.  LIDO  multi-joint  testing  unit. 

On  each  testing  day  subjects  reported  to  the  laboratory  in  a  fasting  condition  (food  was 
restricted  3  hours  prior  to  evaluation).  For  all  upper  extremity  tests,  the  subject  was 
positioned  so  that  the  axis  of  the  joint  was  directly  in  line  with  the  axis  of  the 
dynamometer  goniometer.  All  measures  were  taken  without  gravity  compensation. 
Maximum  isometric  contraction  (MIC)  measures  were  taken  first.  Dynamometer 
attachments  were  selected  and  placed  to  isolate  the  joint;  the  subject  was  positioned 
on  the  instrument  and  maximally  stabilized;  and  then  the  joint  was  positioned  at  a 
specific  angle.  The  subject  was  instructed  each  time  to  give  maximum  efforts.  The 
subject  performed  three  submaximal  contractions  followed  by  one  MIC.  The  subject 
was  given  3  minutes  of  rest  and  this  procedure  was  repeated  in  the  opposing 
direction.  Next,  the  subject  recovered  for  5  minutes  and  then  performed  the  isokinetic 
testing.  The  subject  was  instructed  to  give  maximum  efforts  for  each  repetition  and  to 
move  through  the  entire  range  of  motion  as  rapidly  and  as  forcefully  as  possible.  The 
4  joint  velocity  settings  (60,  120,  180,  and  240  deg/sec)  were  randomly  assigned.  The 
subject  performed  three  submaximal  contractions  at  the  designated  velocity  followed 
by  five  maximum  contractions.  The  subject  was  given  3  minutes  of  recovery  between 
each  velocity  setting. 


Shoulder  Flexion  and  Extension 

Subject  was  in  a  supinated  position  and  was  stabilized  with  velcro  straps  at  the  waist 
and  across  the  clavicle  just  proximal  to  the  shoulder  to  keep  the  right  shoulder  firmly  in 
contact  with  the  plinth  (Figure  2).  The  angle  for  isometric  testing  was  90  degrees.  A  cuff 
attached  to  the  dynamometer  arm  was  placed  on  the  upper  arm  just  proximal  to  the 
elbow  joint.  Therefore,  the  shoulder  joint  was  isolated  and  force  was  applied  at  the 
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point  of  the  cuff  attachment.  The  range  of  motion  where  shoulder  flexion  and  extension 
torque  was  measured  was  between  20  and  1 80  degrees  of  shoulder  flexion. 

Shoulder  Abduction  and  Adduction 

Subject  was  lying  on  the  side  and  was  stabilized  with  velcro  straps  at  the  pelvis  and 
across  the  upper  chest  at  the  axilla  line  to  keep  the  subject's  chest  wall  firmly  in 
contact  with  the  plinth  (Figure  3).  The  angle  for  isometric  testing  was  90  degrees.  A 
cuff  attached  to  the  dynamometer  arm  was  placed  on  the  upper  arm  just  proximal  to 
the  elbow  joint.  Therefore,  the  shoulder  joint  was  isolated  and  force  was  applied  at  the 
point  of  the  cuff  attachment.  The  range  of  motion  where  shoulder  abduction  and 
adduction  torque  was  measured  was  between  15  and  145  degrees  of  shoulder 
abduction. 


Figure  2.  Shoulder  flexion  and  extension. 
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Figure  3.  Shoulder  abduction  and  adduction. 


Shoulder  Internal  (Medial)  and  External  (Lateral)  Rotation 

Subject  was  in  a  supinated  position  with  the  shoulder  placed  at  90  degrees  of 
abduction.  The  subject  was  stabilized  with  velcro  straps  at  the  waist  and  across  the 
mid-upper  arm  to  keep  the  right-upper  arm  and  elbow  firmly  in  contact  with  the  plinth 
(Figure  4).  The  angle  for  isometric  testing  was  0  degrees.  A  cuff  attached  to  the 
dynamometer  arm  was  placed  on  the  forearm  just  distal  to  the  elbow  joint.  The 
shoulder  joint  was  isolated  and  force  was  applied  at  the  point  of  the  cuff  attachment. 
The  range  of  motion  where  shoulder  internal  and  external  rotation  torque  was 
measured  was  between  0  and  70  degrees  of  internal  rotation  and  0  and  60  degrees  of 
external  rotation. 

Elbow  Flexion  and  Extension 

Subject  was  in  a  supinated  position  and  was  stabilized  with  velcro  straps  at  the  waist 
and  across  the  chest  and  mid-upper  arm  to  keep  the  right-upper  arm  and  elbow  firmly 
in  contact  with  the  plinth  (Figure  5).  The  angle  for  isometric  testing  was  60  degrees.  A 
cuff  attached  to  the  dynamometer  arm  was  placed  on  the  forearm  just  proximal  to  the 
wrist  joint.  The  elbow  joint  was  isolated  and  force  was  applied  at  the  point  of  the  cuff 
attachment.  The  range  of  motion  where  elbow  flexion  and  extension  torque  was 
measured  was  between  15  and  135  degrees  of  elbow  flexion. 
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Wrist  Flexion  and  Extension 

Subject  was  in  a  sitting  position  and  was  stabilized  with  velcro  straps  at  the  waist  and 
across  the  chest  to  keep  the  subject's  back  firmly  in  contact  with  the  seat  (Figure  6). 

An  attachment  for  the  forearm  was  secured  to  the  side  of  the  LIDO.  The  subject's  right 
forearm  was  secured  in  a  supinated  position  into  the  device  by  three  small  velcro 
straps:  1)  just  distal  to  the  elbow  joint,  2)  mid-forearm,  and  3)  just  proximal  to  the  wrist 
joint  (Figure  7).  This  configuration  was  designed  to  keep  the  forearm  firmly  in  contact 
with  the  forearm  stabilizing  device.  The  angle  for  isometric  testing  was  0  degrees. 

The  subject  gripped  a  handle  device  attached  to  the  dynamometer  shaft.  The  wrist 
joint  was  isolated  and  force  was  applied  at  the  point  of  the  handle  attachment.  The 
range  of  motion  where  wrist  flexion  and  extension  torque  was  measured  was  between 
0  and  60  degrees  for  wrist  flexion  and  0  and  45  degrees  for  wrist  extension. 


Figure  4.  Shoulder  internal  (medial)  and  external  (lateral  rotation). 
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Figure  5.  Elbow  flexion  and  extension. 


Figure  6.  Wrist  flexion  and  extension. 
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Figure  7.  Closeup  view  of  wrist  flexion  and  extension. 


Wrist  Radial  and  Ulnar  Deviation 

Subject  was  in  a  sitting  position  and  was  stabilized  with  velcro  straps  at  the  waist  and 
across  the  chest  to  keep  the  subject's  back  firmly  in  contact  with  the  seat  (Figure  6). 

An  attachment  for  the  forearm  was  secured  to  the  side  of  the  LIDO.  The  subject's  right 
forearm  was  secured  in  a  supinated  position  into  the  device  by  three  small  velcro 
straps:  1)  just  distal  to  the  elbow  joint,  2)  mid-forearm,  and  3)  just  proximal  to  the  wrist 
joint  (Figure  8).  This  configuration  was  designed  to  keep  the  forearm  firmly  in  contact 
with  the  forearm  stabilizing  device.  The  angle  for  isometric  testing  was  0  degrees 
(neutral).  The  subject  gripped  the  handle  device  attached  to  the  dynamometer  shaft. 
Therefore,  the  wrist  joint  was  isolated  and  force  was  applied  at  the  point  of  the  handle 
attachment.  The  range  of  motion  where  wrist  radial  and  ulnar  deviation  torque  was 
measured  was  between  0  and  35  degrees  for  wrist  radial  deviation  and  0  and 
35  degrees  for  wrist  ulnar  deviation. 

Wrist  (Forearm)  Pronation  and  Supination 

Subject  was  in  a  sitting  position  and  was  stabilized  with  velcro  straps  at  the  waist  and 
across  the  chest  to  keep  the  subject's  back  firmly  in  contact  with  the  seat  (Figure  6). 

An  attachment  for  the  forearm  was  secured  to  the  side  of  the  LIDO.  The  subject's  right 
forearm  was  in  a  neutral  position  and  2  small  velcro  straps  were  placed:  1)  just  distal 
to  the  elbow  joint  and  2)  mid-forearm  (Figure  9).  This  configuration  was  designed  to 
keep  the  proximal  section  of  the  forearm  firmly  in  contact  with  the  forearm  stabilizing 
device  and  to  allow  pronation  and  supination  to  occur.  The  angle  for  isometric  testing 
was  0  degrees  (neutral).  The  subject  gripped  the  handle  device  attached  to  the 
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dynamometer  shaft.  Therefore,  the  forearm  was  isolated  and  force  was  applied  at  the 
point  of  the  handle  attachment.  The  range  of  motion  where  wrist  pronation  and 
supination  torque  was  measured  was  between  0  and  60  degrees  for  wrist  pronation 
and  from  0  to  60  degrees  for  wrist  supination. 


Figure  8.  Wrist  radial  and  ulnar  deviation. 


Figure  9.  Wrist  (forearm)  pronation  and  supination. 


Hip  Flexion  and  Extension 

Subject  was  in  a  supinated  position  and  was  stabilized  with  velcro  straps  at  the  chest 
and  waist  and  across  the  left-upper  thigh  to  keep  the  back  and  pelvis  firmly  in  contact 
with  the  plinth  (Figure  10).  The  angle  for  isometric  testing  was  90  degrees.  A  large 
cuff  attached  to  the  dynamometer  arm  was  placed  on  the  thigh  just  proximal  to  the 
knee  joint.  Therefore,  the  hip  joint  was  isolated  and  force  was  applied  at  the  point  of 
the  cuff  attachment.  The  range  of  motion  where  hip  flexion  and  extension  torque  was 
measured  was  between  0  and  110  degrees  of  hip  flexion. 

Hip  Abduction  and  Adduction 

Subject  was  lying  on  the  side  and  was  stabilized  with  velcro  straps  at  the  chest  and 
waist  and  across  the  left-upper  thigh  to  keep  the  subject's  chest  wall  and  pelvis  firmly 
in  contact  with  the  plinth  (Figure  1 1).  The  angle  for  isometric  testing  was  0  degrees.  A 
large  cuff  attached  to  the  dynamometer  arm  was  placed  on  the  thigh  just  proximal  to 
the  knee  joint.  Therefore,  the  hip  joint  was  isolated  and  force  was  applied  at  the  point 
of  the  cuff  attachment.  The  range  of  motion  where  hip  abduction  and  adduction  torque 
was  measured  was  between  0  and  60  degrees  of  hip  abduction. 

Hip  Internal  (Medial)  and  External  (Lateral)  Rotation 

Subject  sat  on  an  elevated  table  and  was  stabilized  with  velcro  straps  at  the  waist  and 
across  the  left-upper  thigh  to  keep  the  subject's  pelvis  firmly  in  contact  with  the 
exterior  table.  The  angle  for  isometric  testing  was  0  degrees  (neutral).  A  cuff  attached 
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to  the  dynamometer  arm  was  placed  on  the  lower  leg  just  distal  to  the  knee  joint. 
Therefore,  the  hip  joint  was  isolated  and  force  was  applied  at  the  point  of  the  cuff 
attachment.  The  range  of  motion  where  hip  internal  and  external  rotation  torque  was 
measured  was  between  0  and  1 5  degrees  of  internal  rotation  and  from  0  to 
35  degrees  of  external  rotation. 


Figure  10.  Hip  flexion  and  extension. 


Figure  11.  Hip  abduction  and  adduction. 
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Knee  Flexion  and  Extension 

Subject  was  in  a  sitting  position  and  was  stabilized  with  velcro  straps  at  the  waist  and 
across  the  chest  to  keep  the  subject's  back  firmly  in  contact  with  the  seat  (Figure  12). 
The  subject's  right  thigh  was  stabilized  with  an  attachment  secured  to  the  side  of  the 
LIDO.  This  configuration  was  designed  to  keep  the  thigh  firmly  secured  between  the 
stabilizing  device  and  the  seat.  The  angle  for  isometric  testing  was  60  degrees.  A  cuff 
attached  to  the  dynamometer  arm  was  placed  on  the  lower  leg  just  proximal  to  the 
ankle  joint.  Therefore,  the  knee  joint  was  isolated  and  force  was  applied  at  the  point 
of  the  cuff  attachment.  The  range  of  motion  where  knee  flexion  and  extension  torque 
was  measured  was  between  5  and  100  degrees  of  knee  flexion. 

Ankle  Plantarflexion  and  Dorsiflexion 

Subject  was  in  a  supinated  position  with  the  knee  joint  fully  extended.  The  subject  was 
stabilized  with  velcro  straps  at  the  chest  and  waist  to  keep  the  subject's  back  and 
pelvis  firmly  in  contact  with  the  plinth  (Figure  13).  The  subject's  right  knee  was 
stabilized  with  an  attachment  secured  to  the  side  of  the  LIDO.  This  configuration  was 
designed  to  keep  the  knee  joint  firmly  secured  between  the  stabilizing  device  and  the 
plinth.  The  angle  for  isometric  testing  was  0  degrees  (neutral).  The  subject's  right  foot 
was  secured  by  three  small  velcro  straps  and  placed  into  a  foot  plate  device 
(Figure  14).  The  foot  plate  device  was  attached  to  the  dynamometer  shaft.  Therefore, 
the  ankle  joint  was  isolated  and  force  was  applied  at  the  foot  plate  attachment.  The 
range  of  motion  where  ankle  plantarflexion  and  dorsiflexion  torque  was  measured 
was  between  0  and  30  degrees  for  ankle  plantarflexion  and  from  0  to  10  degrees  for 
ankle  dorsiflexion. 


Figure  12.  Knee  flexion  and  extension. 
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Figure  13.  Ankle  plantarflexion  and  dorsiflexion. 


Figure  14.  Closeup  view  of  ankle  plantarflexion  and  dorsiflexion. 


Simulated  Ratchet  Pushing  and  Pulling  Maneuver  Torque 

The  axis  of  the  dynamometer  goniometer  was  positioned  in  line  with  the  subject's 
greater  trochanter.  For  all  simulated  ratchet  pushing  and  pulling  maneuvers,  the 
subject  was  in  a  sitting  position.  The  right  upper  extremity  was  evaluated.  The  subject 
was  stabilized  with  velcro  straps  at  the  waist  and  across  the  chest  to  keep  the  subject's 
back  firmly  in  contact  with  the  seat  (Figure  15).  The  subject  gripped  a  simulated 
ratchet  device  at  a  height  of  90%  of  the  linear  distance  measured  from  the  subject's 
greater  trochanter  to  the  acromioclavicular  joint  (Figure  16).  The  subject  was 
instructed  to  move  the  ratchet  device  from  a  full  back  position  (shoulder  joint  at 
10.9  +/-  0.8  degrees;  elbow  joint  at  1 12.0  +/-  6.0  degrees)  to  a  full  forward  position 
(shoulder  joint  at  66.0  +/-  6.8  degrees;  elbow  joint  at  0  degrees)  [ratchet  push 
maneuver],  and  then  to  return  to  the  starting  full  back  position  [ratchet  pull  maneuver]. 
The  subject  was  told  to  give  maximum  efforts  for  each  repetition  and  to  move  through 
the  entire  range  as  rapidly  and  forcefully  as  possible,  but  to  not  allow  the  scapula 
(upper  back)  to  lose  contact  with  the  seat.  Individual  subject  measures  of  ratchet  grip 
height  and  shoulder  and  elbow  joint  angle  of  excursion  were  recorded.  For  all 
ratcheting  conditions,  the  subject  performed  three  submaximal  contractions  followed 
by  five  maximum  contractions.  The  subject  was  given  5  minutes  of  recovery  between 
each  of  the  3  ratcheting  conditions.  The  first  ratcheting  condition  consisted  of  two 
trials.  The  first  trial  consisted  of  a  constant  resistance  of  25  ft-lbs  (isotonic  load)  for  the 
ratchet  push  maneuver  and  0  ft-lbs  of  resistance  for  the  pull  manuever.  The  subject 
was  given  a  5-minute  recovery  period  and  then  performed  the  second  trial  with  0  ft-lbs 
of  resistance  for  the  ratchet  push  maneuver  and  a  constant  resistance  of  25  ft-lbs  for 
the  pull  maneuver.  The  second  condition  was  the  push  and  pull  ratcheting  done  at  a 
constant  angular  velocity  of  120  deg/sec  (isokinetic).  The  third  condition  was  the  push 
and  pull  ratcheting  done  at  a  constant  angular  velocity  of  240  deg/sec  (isokinetic).  For 
each  repetition,  torque  and  angle  position  data  were  measured. 
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Figure  15.  Front  view  of  simulated  ratchet  pushing  and  pulling  maneuver. 


Figure  16.  Side  view  of  simulated  ratchet  pushing  and  pulling  maneuver. 
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Assessment  of  Anthropometric  Data 

Anthropometric  data  were  gathered  for  height,  weight,  age,  sex,  skin  fold  measures, 
and  dimensional  assessment.  Height  was  measured  in  centimeters  and  weight  in 
kilograms  on  a  physician  scale.  Skin  fold  measures  were  taken  with  a  Lange  skin  fold 
caliper.  Males  were  measured  at  the  abdomen  and  anterior  thigh  and  chest;  and 
females  were  measured  at  the  anterior  thigh,  supraillac  crest,  and  triceps.  Estimation 
of  body  fat  and  lean  body  mass  were  performed  using  equations  specific  for  sex  and 
age  [4].  See  Table  1  for  a  summary  of  this  data. 


TABLE  1 .  Lean  Body  Weight 


SUBJ 

SEX 

AGE 

HT  (cm) 

WT  (kg) 

%Body  Fat 

Lean  Body 

1. 

M 

23 

173 

64.9 

10.4 

58.2 

2. 

M 

25 

178 

76.0 

10.4 

68.1 

3. 

M 

28 

188.3 

84.5 

14.5 

72.2 

4. 

F 

23 

172.3 

7.2 

29.5 

54.4 

5. 

M 

22 

185.0 

88.2 

10.7 

78.8 

6. 

M 

25 

180.0 

86.0 

5.7 

81.1 

7. 

M 

26 

176.3 

95.2 

19.0 

77.1 

8. 

F 

23 

174.0 

60.6 

18.6 

49.3 

9. 

F 

22 

168.0 

59.1 

14.8 

50.4 

10. 

F 

21 

158.5 

46.4 

17.2 

38.4 

11. 

F 

21 

158.0 

51.6 

19.5 

41.5 

12. 

M 

21 

178.5 

80.7 

8.9 

73.5 

13. 

F 

23 

166.0 

55.7 

20.8 

44.1 

14. 

M 

23 

162.5 

68.5 

11.3 

60.8 

The  dimensional  anthropometric  data  were  taken  with  cloth  tape  measures  [1].  The 
time  required  to  make  all  anthropometric  measures  was  0.25  to  0.5  hours  per  subject. 


3.2  Data  Reduction 

A  set  of  streamlined  programs  was  developed  to  process  the  raw  strength  data 
(produced  directly  by  the  LIDO  force  torque  dynamometer)  into  a  compact  polynomial 
coefficient  format.  The  raw  data  were  collected  using  the  LIDOACT  software  executing 
on  an  IBM  PC.  The  files  produced  on  the  PC  were  transferred  to  the  VAX  system  using 
a  data  communication  software  package  (Kermit)  in  binary  mode.  These  data  were 
then  transferred  to  a  UNIX-based  workstation. 

After  being  separated  into  files  by  subject,  velocity,  direction,  and  degree  of  freedom 
for  each  joint,  each  torque  versus  angle  data  file  was  viewed  graphically  and  edited  for 
extraneous  data  points  (Figure  17).  Figure  17  shows  that  the  initial  and  final  portions 
of  the  curves  were  omitted  because  of  the  startup  time  during  which  the  subject  is 
beginning  to  apply  a  maximum  torque.  At  the  end  of  a  motion,  the  subject  was 
anticipating  the  stopping  and  change  of  direction  of  the  LIDO  actuator  arm.  These 
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transition  regions  of  torque  were  inconsistent  and  so  were  not  part  of  our  modeling 
effort. 


Elbow  Flexion 


Raw  Data 


Filtered  Data 


40.0 


Degrees 

Figure  17.  Data  before  and  after  visual  editing. 


75.0 

Degrees 


150.0 


The  following  is  the  flow  of  data  from  raw  files  to  a  torque  function  coefficient  file 
(Figure  18).  This  is  the  actual  UNIX  script  file  used: 

#uncompress  the  large  data  files 
uncompress  *.raw 

"Is"  sh*.raw|  lido  |  tosort  |  sort  |  toffc  >  right_shoulder.ffc 
"Is"  el*.raw|  lido  |  tosort  |  sort  |  toffc  >  right_elbow.ffc 
“Is”  wr*.raw|  lido  |  tosort  |  sort  |  toffc  >  right_wrist.ffc 
#clean  up 
compress  *.raw 
mv  *.xy*  xy 
rm  *.asc 

The  following  is  an  explanation  of  the  above  script  file: 

1 .  The  first  line  of  the  script  file  uncompresses  the  data  files  needed. 

2.  The  "Is"  command  feeds  a  list  of  the  file  names  (one  at  a  time)  to  the  LIDO 
program. 
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3.  The  LIDO  program  converts  each  file  it  receives  from  the  LIDO  format  into  an 
ASCII  format  and  passes  the  data  (one  file  at  a  time)  to  tosort. 

4.  The  tosort  program  computes  the  regression  equations  and  collates  person 
data  and  passes  that  condensed  data  on  to  the  sort  program.  In  addition,  it 
creates  xy  files  of  force  versus  angle  and  stores  them  for  review  later. 

5.  The  sort  program  sorts  the  data  on  each  field  and  passes  that  on  to  the  toffc 
program. 

6.  The  toffc  program  processes  sorted  data  to  produce  files  in  the  torque 
function  coefficient  format  (Figure  18). 

7.  The  original  data  is  then  recompressed  in  the  seventh  line  to  conserve  disk 
space. 

8.  The  xy  files  are  moved  into  a  separate  directory  called  xy  in  the  eighth  line. 

9.  All  the  unnecessary  files  are  cleaned  up  in  the  ninth  line. 

The  following  is  a  quick  description  of  each  of  the  processing  programs: 

lido.c  -  This  program  converted  a  raw  LIDO  generated  data  file  to  a  standard 
output  ASCII  file.  The  LIDO  generated  data  file  (extension  .raw)  had  a  main 
header  area  and  several  subheader  areas.  Lines  in  these  areas  were  terminated 
with  a  line  feed  as  in  a  normal  ASCII  text  file.  The  data  lines  were  terminated  with 
a  carriage  return  and  no  line  feed.  The  data  lines  were  three  integer  values  of 
four  characters  each.  The  first  value  was  the  corrected  torque  at  the  cuff  [5,  6],  the 
second  value  was  the  angle,  and  the  third  was  the  uncorrected  torque  at  the  shaft 
of  the  LIDO  actuator.  The  corrected  torque  values  were  chosen  for  further 
processing. 

tosort.c  -  This  program  converted  the  ASCII  LIDO  data  file  to  torque  function 
coefficients  and  an  associated  subject  name,  joint  name,  direction,  and  velocity. 
The  output  was  used  as  standard  input  for  the  UNIX  sort  utility.  This  sorted  output 
was  then  piped  to  the  toffc.c  program. 

toffc.c  -  This  was  a  program  which  processed  output  from  the  tosort  program  and 
converted  it  into  a  torque  coefficient  file  (Figure  18). 
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right_shoulder  /*  joint  name*/ 
x  /*  axis  */ 

abduction  /*  direction*/ 

4  /*number  of  velocities*/ 

/*  1  velocity  ,  3  polynomial  coefficients,  note:  Y  =  A  +  Bx  +  Cx**2  where  A,B>C  are  the  coefficients*/ 


60.000000 

120.000000 

180.000000 

240.000000 

adduction 

4 

60.000000 

120.000000 

180.000000 

240.000000 

y 

extension 

4 

60.000000 

120.000000 

180.000000 

240.000000 

flexion 

4 

60.000000 

120.000000 

180.000000 

240.000000 

z 

lateral 

4 

60.000000 

120.000000 

180.000000 

240.000000 

medial 

4 

60.000000 

120.000000 

180.000000 

240.000000 


3.477892E+01 
5.058879E+01 
3.441 185E+01 
4.630580E+01 


3.899835E+01 

3.210880E+01 

1.128712E+01 

2.879990E+01 


1.853307E+01 
1 .598296E+01 
1.51 601 0E+01 
9.1 0331 1E+00 


6.689232E+01 

4.769003E+01 

5.670030E+01 

4.729467E+01 


1.650289E+01 
1.434598E+01 
1 .8731 10E+01 
1.489092E+01 


2.395702E+01 
2.1 24104E+01 
2.144591E+01 
1.980874E+01 


2.220639E-01 

-1.974063E-01 

-2.165272E-02 

-2.306546E-01 


-2.004480E-02 
6.244854E-02 
4.01 3751 E-01 
-4.25531 2E-02 


1.233500E-01 
2.713108E-01 
2.109192E-01 
1. 32561 5E-01 


-4.941 486E-01 
-2.906334E-01 
-4.014365E-01 
-3.183639E-01 


-4.517265E-02 

-2.340409E-01 

-1.003202E-01 

-4.907729E-02 


1 .1 28039E-01 
1.564447E-01 
1 .82661 7E-01 
1 .33621 7E-01 


-2.821 324E-03 
-2.446309E-04 
-1.170670E-03 
6.999267E-05 


-1.001798E-04 
-3.269533E-04 
-1.763708E-03 
6.87231 5E-04 


-1 .1 50867E-04 
-9.975418E-04 
-8.056303E-04 
-4.60291 2E-05 


1 .258664E-03 
8.198689E-04 
1.149240E-03 
8.940246E-04 


-5.588200E-04 
-2.806892E-03 
-1 .389751 E-03 
4.31 1464E-04 


-1.947907E-04 
-3.1 86225E-03 
-1.797236E-03 
-8.705539E-04 


Figure  18.  Example  of  the  coefficient  file  format. 


21 


3.3  Regression  Model  Development 


Creating  the  regression  model  involved  processing  the  data  from  each  joint,  direction 
and  velocity,  into  separate  files.  There  were  three  phases  in  the  processing  of  these 
files  which  allow  the  prediction  of  individual  torque  curves. 

1 .  Creating  lean  body  mass  to  average  torque  relationship  tables. 

2.  Generating  a  population  shape  coefficient  table. 

3.  Linking  the  population  shape  table  to  the  lean  body  mass  average  torque 
table. 

Lean  Body  Mass  to  Average  Torque  Table  Generation 

A  relationship  exists  between  lean  body  mass  and  average  torque.  Figure  19  shows 
representative  examples  of  this  linear  relationship  for  the  shoulder,  wrist,  and  elbow. 
Correlation  coefficient  values  range  from  0.80  to  0.95.  The  lean  body  mass  for  each 
of  the  subjects  was  calculated  from  the  measured  weight  and  %  body  fat: 

lean  body  mass  =  weight  *  (1  -  %  body  fat) 

The  average  torque  was  then  computed  for  the  entire  joint  range  for  each  velocity  of 
each  joint.  From  these  data,  linear  regression  coefficients  were  obtained  and 
organized  into  coefficient  tables.  These  tables  can  be  used  to  calculate  the  average 
torque  for  a  particular  joint  and  velocity  when  the  lean  body  mass  is  known. 

Obtaining  Population  Curves 

In  our  model  an  assumption  was  made  that  for  an  isolated  joint  motion,  individuals  use 
the  same  muscle  groups.  These  muscle  groups  basically  have  the  same  shape, 
orientation,  and  points  of  attachment  and  differ  in  the  magnitude  of  force  exerted. 
Hence,  when  the  torque-position  curves  are  normalized  and  plotted,  a  general  trend 
can  be  seen  for  a  particular  joint.  Figure  20  represents  the  normalized  data  for  all  the 
subjects  for  the  elbow  and  shoulder  flexion.  These  data  show  a  definite  trend.  To 
take  advantage  of  this  relationship,  the  data  were  processed  in  the  following  way: 

1 .  We  normalized  each  individual  data  set  (isolated  joint  for  each  axis, 
direction,  and  velocity)  with  respect  to  the  average  torque  for  that  data  set. 

2.  For  all  subjects,  the  normalized  data  for  an  isolated  joint  (axis,  direction, 
and  velocity)  were  combined  into  one  file.  A  total  of  7  axes,  2  directions, 
and  4  velocities  resulted  in  56  composite  files  of  normalized  subject  data. 
(Figure  20  shows  a  sample  case.) 

3.  A  regression  for  each  of  these  combined  files  was  then  computed. 

4.  The  second  order  polynomials,  representing  the  shapes  of  the  joint  motions, 
were  then  organized  into  coefficient  tables  (Figure  18).  These  tables  allow 
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the  computation  of  a  normalized  torque-position  curve  for  any  joint,  axis, 
direction,  and  velocity. 

Generation  of  an  Individual's  Torque-Position  Curve 

To  generate  a  torque-position  curve  from  lean  body  mass  information  it  is  necessary  to 
multiply  the  normalized  torque-position  curve  coefficients  (established  in  the  table  from 
step  2)  by  the  average  torque.  The  average  torque  is  computed  from  the  data  in  the 
lean-body-mass  to  average-torque-coefficient  table  established  in  step  1  above. 

These  two  tables  provide  the  data  to  compute  a  polynomial  which  represents  an 
individual's  torque-position  curve  for  that  joint,  axis,  direction,  and  velocity. 

For  all  of  our  data,  torque  was  measured  for  four  separate  velocities.  No  attempt  was 
made  to  find  a  correlation  between  these  velocities.  Velocity  dependence  was  not 
consistent  among  all  joints.  Since  dynamic  strength  is  important  in  most  joint  motions, 
we  feel  that  this  area  requires  more  attention. 
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Figure  19.  Lean  body  mass  (kg)  versus  mean  torque  (ft-lbs). 
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Wrist  Flexion 


Wrist  Radial  Deviation 


w(1-*bf) 


w(  1  -%bf) 


Wrist  Supination 
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Figure  19.  Lean  body  mass  (kg)  versus  mean  torque  (ft-lbs)  (continued). 
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Elbow  Flexion:  velocity  =  120  deg/sec 


Normalized  Data:  1 4 Subject®  Polynomial  Fit  to  Normalized  Data 


Degrees  Degrees 


Shoulder  Flexi 
Normalized  Data:  1 4  Subjects 

2.6 


Torque 
1 .3 


0.0 

0.0  100.0  200.0 


on:  velocity  =  120  deg/sec 

Polynomial  Fit  to  Normalized  Data 


Degrees 


Degrees 


Figure  20.  Polynomial  regression  fit  to  normalized  data. 
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4.0  RESULTS  AND  DISCUSSION 


Fourteen  subjects  were  measured  for  isolated  joint  torques  in  each  axis,  direction,  and 
for  four  different  velocities.  This  study  shows  how  these  time-consuming  measure¬ 
ments  of  dynamic  isolated  joint  torques  can  be  predicted.  A  strong  correlation 
between  torque  produced  by  an  isolated  joint  and  the  individual's  lean  body  mass  was 
found.  This  relationship  allows  prediction  of  the  torque-position  curves  for  each 
isolated  joint  by  a  simple  index,  the  lean  body  mass. 

To  verify  our  regression  model  for  predicting  isolated  joint  torque-position  curves,  we 
plotted  predicted  versus  measured  torque-position  curves.  Figure  21  is  an  example  of 
a  prediction  of  an  isolated  joint  torque-position  curve  (from  lean  body  mass  information 
only)  plotted  with  measured  data.  The  general  shape  of  the  curves  of  the  data  are 
reproduced. 


Velocity  =  120  deg/sec 


Shoulder  Flexion 


Elbow  Flexion 


Figure  21.  Isolated  joint  curve:  predicted  versus  measured. 


To  get  an  estimate  of  the  deviation  from  the  measured  data,  a  comparison  computation 
was  done.  For  each  of  the  torque-position  data  sets  collected,  a  corresponding  array 
of  torque  values  was  calculated  from  the  predicted  polynomial  coefficients  for  that 
particular  data  set.  A  torque  difference  vector  (a  difference  of  the  actual  measured 
torque  value  at  that  angle  and  the  computed  value  at  that  point  from  the  predicted 
polynomial  coefficients)  was  then  created.  A  percentage  absolute  value  (relative  to 
the  maximum)  of  the  difference  array  was  calculated  and  plotted.  This  same  analysis 
was  done  on  the  regression  coefficients  computed  from  that  particular  data  set.  The 
regression  coefficients  represent  the  ideal  curve  through  the  collected  data  values. 

Figures  22,  23,  and  24  represent  a  comparison  of  the  error  of  the  measured  data  fit 
regression  coefficients  versus  the  error  from  the  lean  body  mass  predicted  coefficients 
for  one  representative  subject.  The  error  of  the  regression  coefficients  reflects  the 
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deviation  in  the  fit  to  the  measured  data.  In  these  figures,  the  average  error  of  the  fit  to 
the  measured  data  is  6%.  For  prediction  of  second  order  regression  coefficients  from 
lean  body  mass  information,  this  is  the  best  that  one  can  hope  to  achieve.  As 
expected,  the  predicted  coefficients  from  lean  body  mass  have  a  greater  deviation 
than  the  regression  coefficients  calculated  from  the  actual  data.  For  the  shoulder 
flexion/extension,  elbow  flexion/extension,  and  wrist  flexion/extension  the  deviations 
are  9.8%  for  this  individual.  Figure  25  represents  the  same  comparison  for  all  joints 
for  all  individuals.  When  compared  across  the  entire  subject  pool,  the  shoulder  and 
elbow  joint  prediction  equation  represents  about  a  12%  absolute  value  deviation. 

The  wrist  joint  has  a  greater  deviation.  The  wrist  data  have,  on  the  average,  one-tenth 
the  range  of  the  shoulder  and  elbow  joints.  As  a  result  of  this  small  range,  a 
comparatively  small  deviation  translates  into  a  large  percentage  difference.  The  wrist 
data  need  further  consideration. 

If  only  lean  body  mass  data  on  a  particular  individual  are  available,  our  results  indicate 
that  the  prediction  equations  can  yield  results  on  the  order  of  12%  deviation  from  the 
actual  data  for  all  axes  and  directions  of  the  elbow  and  shoulder  joints.  The  results 
here  could  also  be  used  to  fill  in  missing  or  partial  strength  information  for  an 
individual.  For  instance,  if  only  one  axis  of  a  particular  joint  was  measured,  the  other 
strength  information  of  the  other  axes  and  directions  could  be  extrapolated  from  our 
equations  taking  into  account  the  known  information.  Other  correlates  to  dynamic 
strength  prediction  may  exist.  An  area  of  future  research  is  finding  multiple  correlates 
of  strength. 
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Shoulder  Flexion/Extension 

I  Predicted  from  lean  body  mass 
H  Regression  of  measured  data 


60  120  180  240  60  120  180  240 


Extension - 1 - Flexion 

Velocity  (deg  /sec) 


Figure  22.  One  subject's  shoulder  joint  absolute  value  error  of  the  predicted  curve  as 
compared  to  the  error  of  the  regression  curve  generated  through  the  collected  data. 
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Elbow  Flexion  /Extension 


20  n 


■  Predicted  from  lean  body  mass  only 
I  Regression  of  measured  data 


- Extension - 1 - Flexion 

Velocity  (deg/sec) 


Figure  23.  One  subject's  elbow  joint  absolute  value  error  of  the  predicted  curve  as 
compared  to  the  error  (variation)  of  the  regression  curve  generated  through  the 

collected  data. 
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I  Predicted  from  lean  body  mass 
H  Regression  of  measured  data 
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Figure  24.  One  subject's  wrist  joint  absolute  value  error  of  the  predicted  curve  as 
compared  to  the  error  (variation)  of  the  regression  curve  generated  through  the 

collected  data. 
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Figure  25.  Entire  subject  pool.  This  is  the  absolute  value  error  of  the  predicted  curve 
as  compared  to  the  error  (variation)  of  the  regression  curve  for  all  subjects  for  all  axes 

of  the  joints. 


Our  model  was  based  on  measurements  of  14  subjects.  These  subjects  were  all 
healthy  young  adults.  Extending  our  lean  body  mass-torque  regression  model  beyond 
this  scope  would  require  additional  validation. 

Our  results  include  several  useful  tables: 

1 .  Tables  of  second  order  polynomial  coefficients  were  generated  from  the 
measurements  for  each  joint,  axis,  direction,  and  velocity.  These  tables 
provide  a  convenient  and  compact  storage  and  retrieval  mechanism  to 
access  our  entire  measured  subject  pool  data  of  dynamic  isolated  joint 
torques  (Appendix  A). 

2.  Tables  of  first  order  polynomial  coefficients  relating  lean  body  mass  of  an 
individual  to  the  average  torque  for  an  isolated  joint  motion  were  generated. 
They  allow  the  calculation  of  an  average  torque  for  a  particular  joint, 
direction,  and  velocity  from  an  individual's  lean  body  mass  (Appendix  B). 
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3.  Also  available  are  second  order  coefficient  tables  characterizing  the  torque- 
position  curves  that  have  been  normalized  across  our  subject  pool. 

These  curves  represent  the  general  (normalized)  trends  in  a  torque- 
position  curve  and  may  be  useful  in  studying  torque  as  a  function  of  joint 
angle  in  a  population  (Appendix  C). 

To  make  efficient  and  practical  use  of  all  the  coefficients,  variables,  and  relationships 
found  in  this  report,  we  have  encapsulated  the  results  of  this  study  in  a  tool  kit  of 
software  routines  executing  on  a  variety  of  platforms  such  as  UNIX,  DOS,  and 
Macintosh  machines  (Appendix  D).  This  code  contains  our  isolated  joint  torque  model 
(with  tables  of  torque  coefficients)  that  allows  the  prediction  of  dynamic  isolated  joint 
torques  from  an  individual’s  lean  body  mass. 

Although  our  report  has  focused  on  isolated  joint  prediction,  we  did  measure  and 
correlate  a  multi-joint  task  (the  ratchet  wrench  push-pull)  to  lean  body  mass.  Figure  26 
indicates  that  multi-joint  tasks  may  also  be  related  to  lean  body  mass.  Percentage 
body  fat  is  a  good  predictor  of  torque.  There  is  a  strong  correlation  (r  >  0.92)  between 
torque  production  capability  and  the  lean  body  mass.  Once  a  representative  sample 
of  a  population  has  been  measured  for  a  particular  composite  motion,  joint  strength 
and  prediction  of  torque  capability  of  a  particular  individual  may  be  extrapolated  by 
only  two  measures:  percentage  body  fat  and  weight.  Research  continues  in  this  area. 

In  addition,  this  ratchet  data  as  well  as  all  the  isolated  joint  data  are  the  basis  of  a 
graphically  based  human  strength  model  being  developed  at  the  NASA-JSC  Man- 
Systems  Division.  This  model  calculates  end  effector  strength  based  on  any  arbitrary 
motions  using  this  dynamic  isolated  joint  information  and  a  vector  sum  algorithm. 
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Lean  body  mass  vs.  torque  for  ratcheting  task  (push) 


Lean  body  mass  vs.  torque  for  ratcheting  task  (pull) 


Figure  26.  Lean  body  mass  versus  mean  torque  for  ratchet  wrench  maneuver. 
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5.0  FUTURE  DIRECTIONS 

1 .  Extend  isolated  joint  measurements  to  all  joints. 

2.  Extend  the  subject  pool  to  include  a  more  diverse  group. 

3.  Establish  a  database  of  dynamic  isolated  joint  torques  for  general  use. 

4.  Examine  velocity  dependence  in  more  detail. 
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7.0  APPENDICES 


Note:  All  data  and  code  presented  in  this  paper  are  available  in 
electronic  form.  Contact  us. 


APPENDIX  A  -  Polynomial  coefficients  of  angle  versus  torque  collected  data  for  one 
subject. 
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APPENDIX  B  -  Coefficient  table  relating  lean  body  mass  to  average  torque  for  a  particular  velocity. 
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radial 

4 

60.0000 

120.0000 

180.0000 

240.0000 

ulnar 

4 

60.0000 

120.0000 

180.0000 

240.0000 

y 

extension 

4 

60.0000 

120.0000 

180.0000 

240.0000 

flexion 

4 

60.0000 

120.0000 

180.0000 

240.0000 

z 

pronation 

4 

60.0000 

120.0000 

180.0000 

240.0000 

supination 

4 

60.0000 

120.0000 

180.0000 

240.0000 


-1.268475E+00 
-1.005637E+00 
-1.270290E+00 
-1. 77091 2E+00 


-2.900480E+00 

-3.978939E+00 

-2.994360E+00 

-3.975077E+00 


-1. 127101  E+00 
-7.367646E-01 
-1.423397E+00 
-1.275451  E+00 


4.74761 6E+00 
-4.260902E+00 
-2.04081 7E+00 
-3.654124E+00 


-6.592008E-01 

-1.366900E+00 

-1.853492E+00 

-2.182193E+00 


-2.369235E+00 
-2.446581  E+00 
-2.280463E+00 
-2.641 539E+00 


1 .357690E-01 
1.208873E-01 
1.214390E-01 
1 .26391 8E-01 


1.932900E-01 

2.042846E-01 

1.823847E-01 

1.952298E-01 


8.41 6442E-02 
7.451715E-02 
8.59421 3E-02 
8.026770E-02 


2.070452E-01 

1.903496E-01 

1.303384E-01 

1.715256E-01 


7.739327E-02 

9.181561E-02 

1.013228E-01 

1.070073E-01 


1.065645E-01 
1.088279E-01 
1.070542E-01 
1.11 9357E-01 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 
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Appendix  C  -  Normalized  torque  coefficients  for  entire  subject  pool. 
right_elbow 

y 

extension 


4 

60.0000 
120.0000 
180.0000 
240.0000 
flexion 
4 

60.0000 
120.0000 
180.0000 
240.0000 
right_shoulder 
x 

abduction 
4 

60.0000 
120.0000 
180.0000 
240.0000 
adduction 
4 

60.0000 
120.0000 
180.0000 
240.0000 

y 

extension 
4 

60.0000 
120.0000 
180.0000 
240.0000 
flexion 
4 

60.0000 
120.0000 
180.0000 
240.0000 
z 

lateral 
4 

60.0000 
120.0000 
180.0000 
240.0000 
medial 
4 

60.0000 
120.0000 
180.0000 
240.0000 
right_wrist 
x 

radial 


2.382654E-01 

3.998269E-01 

6.186674E-01 

6.797335E-01 


9.436899E-01 
9.559348E-01 
1.024970E+00 
1.181 190E+00 


1. 63821 7E+00 
1.548579E+00 
1.754382E+00 
1.667307E+00 


8.051 392E-01 
8.1 1 4389E-01 
8.734385E-01 
7.575982E-01 


6.296859E-01 

4.782509E-01 

7.084687E-01 

6.559435E-01 


1.873006E+00 
1  857643E+00 
1 .801 71 9E+00 
1 .71 0937E+00 


1 .061 949E+00 
1.041042E+00 
1.007360E+00 
9.669892E-01 


1.022044E+00 
1 .041 949E+00 
1.029437E+00 
1.021 175E+00 


1.609246E-02 
1.352749E-02 
7.465366E-03 
5.1 23406E-03 


1.385255E-03 
2.457487E-03 
7.60521 7E-04 
-2.677844E-03 


-9.23421 3E-03 
-6.3281 32E-03 
-1 .1 59083E-02 
-9.9094 23E-03 


4.1 59272E-03 
4.182187E-03 
1 .71 7887E-03 
4.31 9388E-03 


7.588339E-03 
9.29431 4E-03 
4.1 27454E-03 
3.129994E-03 


-1 .086977E-02 
-1.102275E-02 
-1.150634E-02 
-9.7885 73E-03 


-5.109649E-03 
-4.551 276E-03 
-4.989951  E-03 
-4.986339E-03 


4.982772E-03 

4.882844E-03 

4.798532E-03 

3.547545E-03 


-7.057701  E-05 
-6.401 985E-05 
-2.948732E-05 
-1.389026E-05 


-7.300826E-06 
-1 .9901 86E-05 
-1 .1 20966E-05 
3.1 23384E-06 


1.584508E-05 
-3.656676E-06 
2.1 8831 2E-05 
1.595937E-05 


-1 .9231 75E-05 
-2.022592E-05 
-2.843898E-06 
-1.588424E-05 


-3.253375E-05 
-3.527058E-05 
-1 .1471 1 6E-05 
-1 .1 57555E-06 


2.096022E-05 

2.299863E-05 

2.964644E-05 

2.236658E-05 


-1  644733E-05 
-1.793456E-05 
9.01 7851 E-06 
2.867537E-05 


-3.170801  E-05 
-5.522752E-05 
-4.722543E-05 
-3.615252E-05 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 
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4 

60.0000 

120.0000 

180.0000 

240.0000 

ulnar 

4 

60.0000 

120.0000 

180.0000 

240.0000 

y 

extension 

4 

60.0000 

120.0000 

180.0000 

240.0000 

flexion 

4 

60.0000 

120.0000 

180.0000 

240.0000 

z 

pronation 

4 

60.0000 

120.0000 

180.0000 

240.0000 

supination 

4 

60.0000 

120.0000 

180.0000 

240.0000 


1.181802E+00 

1.107246E+00 

1.036549E+00 

1.039976E+00 


-1.675913E-02 
-1 .877522E-02 
-1 .984581 E-02 
-2.084383E-02 


1.066880E+00 
1.087365E+00 
1.1 14945E+00 
1 .07331 1E+00 


9.806204E-01 

9.499655E-01 

9.177343E-01 

9.216082E-01 


1.084562E+00 
1.109554E+00 
1. 091251 E+00 
1.088930E+00 


1.01 0391  E+00 
9.950561  E-01 
1.009005E+00 
9.778854E-01 


1.024255E+00 
9.929054E-01 
1.009055E+00 
9.361 496E-01 


9.789578E-03 
1 .1 14404E-02 
1.416128E-02 
1.608139E-02 


8.9601 70E-03 
9.603754E-03 
8.428668E-03 
9.430406E-03 


4.1 2461 7E-04 
5.107403E-04 
-5.016060E-05 
-9.37731 2E-04 


5.380042E-03 
5.733134E-03 
5.072503E-03 
6.1 241 77E-03 


-7.277947E-03 
-7.556385E-03 
-7.271 954E-03 
-6.8241 19E-03 


-6.016592E-04 

-4.386299E-04 

-3.919646E-04 

-3.424259E-04 


-3.1 74051 E-04 
-3.558657E-04 
-3.47691 9E-04 
-3.90451 5E-04 


-5.493223E-05 
-4.643842E-05 
-3.641 872E-06 
-2.801 543E-05 


-1.078120E-04 
-1 .394872E-04 
-1.093535E-04 
-1 .0801 13E-04 


3.692308E-06 
7.248831  E-06 
-6.918676E-06 
1.015605E-05 


-3.147846E-05 
-3.043222E-05 
-5.40331 2E-05 
-1 .5931 18E-05 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 


0.000000E+00 

0.000000E+00 

0.000000E+00 

0.000000E+00 
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Appendix  D  -  Source  code  for  utilization  of  torque  coefficient  tables.  This  code  Is  available  compiled 
on  UNIX,  IBM  PC,  and  Macintosh  computers.  Contact  us. 
include  <stdio.h> 

#include  <math.h> 

#include  "lean.h" 

/ . ******* . * . . . . * . ***» . . . — / 

Authors  Abhilash  Pandya,  James  Maida. 

Date:  8/29/91 

Note:  You  are  welcome  to  copy  and  modify  this  set  of  routines  to  fit  your 
own  needs. 

This  is  a  set  of  tool  kit  programs  to  allow  easy  manipulation  of  the 
lean  body  mass  to  torque  relationship  (see  report).  It  creates  coefficient 
files  which  describe  fully  each  joint.  The  coefficient  files  relate  torque 
for  a  particular  joint  to  its  angle  and  velocity.  See  the  discussion  on 
ffc  files  and  the  header  file  which  describes  the  ffc  files. 

The  main  program  is  an  example  of  the  way  the  tool  kit  can  be  utilized. 

It  does  the  following  operations: 

1  ->  Input  a  new  %bf  and  weight 

2  ->  Print  Population  Coefficients 

3  ->  Print  Lean  Body  Mass  coefficients 

4  ->  Print  This  individual’s  coefficients 

5  ->  Write  out  (*.ffc  files)  individuals  coefficients 

6  ->  Print  out  a  torque  summary  for  this  individual 

Units:  weight-  kg.  torque  -  ft-lbs. 

These  routines  include  only  the  standard  C  files  and  should  compile  on 
most  machines  in  its  present  form, 
to  compile  use  : 


cc  lean.c  leanlib.c  4m  -o  lean 
to  execute  use  :  lean  <%body  fat>  <  weight  > 

For  Mac:  Think  C  compiler,  make  sure  to  include  the  ANSI  Library 

in  your  project  file  and  make  sure  that  the  ANSI  library  is  split  into  its 

own  code  segment  ( There  is  a  32  K  code  segment  limit,  see  you  manuals  for  details). 

A************  ***********************  jHHHHHHHHHHHHHHHHHHHHHHHHNHHHHHHHHHHHtr*  *+  ♦  A********** 

//*add  to  this  list  as  more  joints  become  available  and  increment  JC  the  joint  count*/ 
static  char  joint_name[MAXJOINT][80]  =  {"right_shoulder", 

“right_elbow", 

"right_wrisr }; 


static  AxisTable  AxisMapQ  = 

{ 

“ABD  VEL",  "abduction",  "x",  1, 

"ADD  VEL",  "adduction",  "x",-1, 

"FLEX  VEL",  "flexion",  "y",  1, 

"EXT  VEL",  "extension",  Y,-1, 

"EXTROT,  "lateral",  "z",  1, 
"INTROT,  "medial",  "z",-1, 
"PRONATION  VELOCITY",  "pronation", 
"SUPINATION  VELOCITY",  “supination", 
"RADIAL  DEV",  "radial",  V,  1, 

"ULNAR  DEV",  "ulnar",  V, -1, 

0,  0,  0,  0 


N*_ N 

^  I 

^  I 


1, 

-1. 


int  AxisTableSize  =  10; 
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Calculatel ndivid ualCoeffT able  (itable,  ptable,  Itable,  bf,  weight,  count) 

CoeffTable  itableQ; 

CoetfTable  ptableQ; 

CoeffTable  ItableQ; 
float  bf; 
float  weight; 
int  count; 

A*************************************************  ***************************  XT************************************ 

Function: 

Given  the  population  normalized  table  (ptable)  and  the 
lean  body  mass  to  torque  relationship  table  ( Itable ),  this  program  will 
calculate  the  the  itable  ( individual  person's  table)  from  his/her 
body  fat  and  weight. 

The  itable  is  calculated  by  computing  the  average  torque  for  each 
joint  direction  and  velocity  from  the  Itable.  The  Itable  (lean  body  mass 
table)  contains  coefficients  which  describe  the  relationship  between 
the  average  torque  and  the  lean  body  mass.  To  use  this  table,  first 
a  lean  body  mass  is  computed  from  the  weight  and  %  body  fat.  This  value  is 
plugged  into  the  coefficients  looked  up  in  the  Itable.  This  yields  the 
average  torque  for  that  joint  direction  and  velocity.  After  this  torque  is 
computed, the  normalized  torque  coefficients  for  that  joint,  direction  and 
velocity  are  looked  up  in  the  ptable  and  multiplied  by  the  average  torque 
computed  from  the  Itable.  These  coefficients  are  then  stored  into  the  itable 
and  returned. 

Parameters: 

itable  -  individuals  table 

ptable  -  population/  normalized  table 

Itable  -  lean  body  mass  to  average  torque  table 

bf  -  body  fat. 

mass  -  mass  of  the  person. 

***************  ********  **  *******^  **■•*****★★*★*»★★★*★★*  **********  *************  ★»  ************* ********* *********  ***** 

{ 

int  axis; 
int  dir; 

int  eq_count; 
int  i,  j ,  z; 

float  lean_body_mass,  average; 
float  ComputeCoeff  0 

/*lean  body  mass  calculation*/ 
lean_body_mass  =  weight  *  ( 1 .0  -  bf );; 

printf  ( “Coefficients  calculated  for  lean  body  mass  %f ,  lean_body_mass); 

/*for  each  joint*/ 
for  (j  =  0  ;  j  <  count;  j++) 

{ 

/*copy  the  name  and  axis  count  to  the  itable*/ 
strcpy  ( itable[j].joint_name,  ptable[j].joint_name); 
itable[j].axis_count  =  ptable[j].axis_count; 

/*for  each  axis  in  this  joint*/ 

for(  axis  =  0;  axis  <  ptable[j].axis_count;  axis++ ) 

{ 

/*copy  the  axis  number*/ 
itable[j].axis[axis]  =  ptable[j].axis[axis]; 
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/*for  each  direction  in  this  axis*/ 
for(  dir  =  0;  dir  <  2;  dir++ ) 

{ 

/* copy  the  direction  name  and  equation  count  to  the  ptable*/ 
strcpy  ( itable[j].function[axis][dir].direction_name, 
ptable(]].function[axis][dir].direction_name); 
itable[j].function[axis][dir].eq_count  = 
ptable[j].function[axis][dir].eq_count; 
eq_count  =  ptableO].function[axis][dir].eq_count; 

/*for  each  regression  equation  compute  the  individuals  eqns*/ 

for(  i  =  0;  i  <  eq_count;  i++ ) 

{ 

/*copy  the  velocity  information*/ 
itable[j].funcbon[axis][dir].velocity[i]  = 
ptableQ].function[axis][dir].velocity[i]; 

/*compute  the  average  for  torque  for  this  dir,  velocity*/ 

/*from  the  lean  body  mass  information  and  the  Itable  */ 
average  = 

ComputeCoeff  (ltable[j].function[axis][dir].coeffp], 
lean_body_mass); 

/*multiply  the  coefficients  of  the  pop  table  by  average*/ 

/*and  store  in  the  itable  */ 

MultiplyCoeff  (itable[j].function[axis][dir].coeff[i], 
ptableQ].function[axis][dir].coeffIi], 
average); 

}/*endfbri*/ 

}  /*end  for  dir*/ 

}/*  end  for  axis*/ 

}/*end  for  j-joint*/ 


} 

WriteCoeffTableToFile(  table,  ext,  count ) 
CoeffTable  tableQ; 
char  ext[]; 
int  count; 


Function: 

Given  a  table  pointer  a  file  extension,  and  the  number  of  joints,  it  will  write 
the  coefficients  of  the  table  specified  into  files  with  the  names 
of  the  files  taken  form  the  joint  names  in  joint_name  table  with  the  extension 
in  variable  ext. 

This  is  an  example  of  the  file  format: 

right_elbow 

Y 

extension 

4 

60.000000  1.743286e+01  1.1 7741 7e+00  -5.1 6381 8e-03  0.000000e+00 
120.000000  2.869938e+01  9.709965e-01  -4.59531 3e-03  0.000000e+00 
180.000000  4.1 10020e+01  4.959499e-01  -1 ,958944e-03  0.000000e+00 
240.000000  3.991 397e+01  3.008465e-01  -8.156364e-04  0.000000e+00 
extension 
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4 

60.000000  1.743286e+01  1.1 7741 7e+00  -5.1 6381 8e-03  0.000000e+00 
120.000000  2.869938e+01  9.709965e-01  -4.5953136-03  0.000000e+00 
180.000000  4.1 10020e+01  4.959499e-01  -1 .958944e-03  0.000000e+00 
240.000000  3.9913976+01  3.008465e-01  -8.156364e-04  0.000000e+00 

line  1 :  joint  name 
line  2:  axis 
line  3:  direction 

line  4:  number  of  regression  equations, 
line  5:  velocity,  coefficient  1  coefficient  2  ... 
this  repeats  until  the  entire  joint  is  described. 

Parameters: 

table  -  storage  table  for  the  coefficient  values. 

ext  -  file  extension  the  data  to  be  stored  in  ( eg  ffc ,  ntc  or  Ibc ) 

count  -  number  of  joints. 

*»**»**»**»★*★★**»*»»»★**»** ***************** ★ *************** A ****** ★^★★»»»»*»»»**»** »***»»**»**★★*» ************** 

{ 

FILE  *fp; 

char  filename[80j; 
char  line[132]; 

Int  i,  j ; 
int  axis; 

char  axisname[32]; 
char  dirname[80j; 
int  ano; 
int  dir; 

int  eq_count; 

Int  joint; 

/*for  each  joint  in  the  joint  names, 

write  the  file  for  that  joint 

for  the  table  specified.*/ 

for  ( joint  =  0  ;  joint  <  count ;  joint++) 

{ 

/*open  a  file  for  this  joint  using  names  form  the 
joint_name  table 

*/ 

strcpy  ( filename,  joint_nameOoint]); 
sprintf(  filename,  "%s.%s",  filename,  ext  ); 
fp  =  fopen(  filename,  "w" ); 
if(  fp  ==  NULL ) 

{ 

printf(  "failed  open  on  %s\n",  filename ); 
return(0); 

} 

/*  write  out  name  of  joint  */ 
fprintf(  fp,  "%s\n",  table[joint].joint_name  ); 
for  (axis  =  0  ;  axis  <  table[joint].axis_count ;  axis++ ) 

if  ( table[joint].axis[axis]==  0 )  fprintf(fp,  “X\n"); 
if  ( tableDoint].axis[axis]==  1 )  fprintf(fp,  rY\n“); 
if  ( table[joint].axis[axis]==  2 )  fprintf(fp,  "Z\n"); 

/*  direction  name  */ 
fbr(  dir  =  0;  dir  <  2;  dir++ ) 

{ 
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fprintf(fp,  "%s\n“,  table[joint].function[axis][dir].direction_name ); 
/*  number  of  force  equations  for  this  direction  */ 


fprintf(  fp,  “%d\n“, 

table[joint].function[axls][dir].eq_count ); 
eq_count  =  table[joint].function[axis][dir].eq_count; 


} 


for(  i  =  0;  i  <  eq_count;  i++ ) 

{ 

fprintf(  fp,  “%f  %e  %e  %e  %e  \n", 

table[joint].function[axis][dir].velocity[i], 
table[joint].function[axis][dir].coeff[i][0], 
table[jointj.function[axisj[dirj.coeff[ij[lj, 
table[joint].function[axis][dir].coeff[i][2]I 
table[joint].function[axis][dir].coeff[i][3] ); 

} 

} 

}  /*  axis  loop  */ 

fdose  (fp); 

}  /*end  for  each  joint*/ 


TorqueSummary  ( table,  bf,  weight,  count) 

CoeffTable  tableQ; 

float  bf; 

float  weight; 

int  count; 


^*****-*  *********************  **  *  ***********  ********  ************  ***************  *  *****  ***********  *******  *******  ******** 

Function:  Print  out  a  quick-look  table  of  the  average  torque  values 
for  each  of  the  joints  in  the  joint  chain.  This  function  takes  the  lean 
body  mass  table  coefficients  the  %bf,  and  weight  of  an  individual  and  computes 
the  average  torque  for  joint,  direction.  All  velocities  are  averaged  before 
being  printed.  This  means  that  the  dynamic  information  is  not  printed  in  order 
to  keep  it  a  quick  look  table. 

Parameters: 

table  -  the  lean  body  mass  table, 
bf  -%bodyfat(0- 1) 
weight  -  weight  of  the  person, 
count  -  joint  count. 

*******  A****************************************************************************** ***************************** 

{ 

int  axis; 
int  dir; 

int  eq_count; 
int  i,  j ; 

float  average,  avg; 
float  ComputeCoeff  0; 
float  lean_body_mass; 
lean_body_mass  =  weight  *  ( 1 .0  -  bf ); 
printf  (“\n\n“); 

printf  ("  Weight  %f  Body  Fat  %f  Lean  Mass  %f\n",  weight,  bf,  lean_body_mass); 
for  (j  =  0  ;  j  <  count;  j++) 

{ 

printf(  “\n\nJoint  Name:  %s\n",  table[j].joint_name ); 
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for(  axis  =  0;  axis  <  table[j].axis_cx>unt;  axis++ ) 

{ 

for(  dir  =  0;  dir  <  2;  dir++ ) 

{ 

printf(  "Direction:  %s", 

table[j]  .function  [axis]  [dir]  .direction_name ) ; 

eq_count  =  table[j].function[axis][dir].eq_count; 
for(  i  =  0;  i  <  eq_count;  i++ ) 

{ 

/*compute  the  average  for  torque  for  this  dir,  velocity*/ 
avecage  = 

ComputeCoeff  (tableO].function[axis][dir].coeff[i], 
lean_body_mass); 
avg  +=average; 


} 

printf  ( "  Average  Torque:  %An",  avg  /  eq_count); 
avg  =  0.0; 

} 

} 

} 

printf  (" - \nn); 

} 

ReadCoeffTable(  table,  ext,  count ) 

j*  ****  A  **  Hr***  *-**-»**  **  »»  **»»*★*»*  *  ★  ***  ********************  ********  *  A**  *  A*  A*  A  A*  AAA  A  A*  A**-****  A*~A  *  »★*  **★*  ***»  *r*  *  ★**  irk  **rk 

Function: 

Given  a  table  pointer  a  file  extension,  and  the  number  of  joints,  it  will  fill 
the  coefficient  table  specified  with  the  values  in  the  file  according  to 
the  joint  names  in  the  joint_name  table  defined. 

Parameters: 

table  -  storage  table  for  the  coefficient  values. 

ext  -  file  extension  the  data  are  stored  in  (eg  ffc ,  ntc  or  Ibc ) 

count  -  number  of  joints. 

jkkkk  kkkkkkkkkkk+irk*k*kkk++kk1rkk-k+k-k*+1rk+k*1rk^4rtHt++-trk  **********************  Ak-kkk*irkirirkAkAk  A  *********  A  ************** 

CoeffTable  tablet]; 
char  extQ; 
int  count; 

{ 

FILE  *fp; 

char  filename[80]; 
char  line[1 32]; 
int  i,  j ; 
char  *ier; 
int  axis; 

char  axisname[32]; 
char  dirname[80]; 
int  ano; 

int  dir; 

int  eq_count; 
int  joint; 

/*  initialize  force  table  */ 


48 


lnitJointCoeffTable(  table ); 

/* for  each  joint  in  the  joint  names, 
load  the  file  for  that  joint 
into  the  table  specified.*/ 

for  ( joint  =  0  ;  joint  <  count ;  joint++) 

{ 

strcpy  ( filename,  joint_name[joint]); 
sprintf(  filename,  "%s.%s",  filename,  ext  ); 
fp  =  fopen(  filename,  "r" ); 
if(  fp  ==  NULL) 

{ 

printf(  "failed  open  on  %s\n",  filename ); 
return(  0 ); 

} 

/*  read  in  name  of  joint  7 

ier  =  fgets(  line,  132,  fp ); 

if(  ier  ==  NULL )  retum(  0 ); 

sscanf(  line,  "%s\n",  table[joint].joint_name  ); 

/*  read  remaining  data  in  file  7 
axis  =  0; 
while(  1 ) 

{ 

ier  =  fgets(  line,  132,  fp ); 
if(  ier  ==  NULL )  break; 

/*  axis  number  7 

sscanf(  line,  "%s",  axisname  ); 

ano  =  FTGetAxisNumber(  axisnamejO] ); 

if(  ano  ==  -1  ) 

{ 

printf(  "invalid  axis  name  %s\n",  axisname ); 
printf(  "using  default  X  axis\n“ ); 
ano  =  0; 

} 

table[joint].axis[axis]  =  ano; 
for(  j  =  0;  j  <  2;  j++) 

{ 

ier  =  fgets(  line,  132,  fp ); 
if(  ier  ==  NULL )  break; 

/*  direction  name  7 
sscanf(  line,  "%s",  dirname ); 

/*  get  direction  code  0  =  negative,  1  =  positive,  2  =  error  7 
dir  =  GetAxisDirection(  dirname ); 
rf(  dir  ==  2 )  dir  =  0; 

strcpy(  table[joint].function[axis][dir].direction_name, 
dirname ); 

ier  =  fgets(  line,  132,  fp ); 
jf(  ier  ==  NULL )  break; 

/*  number  of  force  equations  for  this  direction  7 
sscanf(  line,  "%d", 

&table[joint].function[axis][dir].eq_count ); 
eq_count  =  table[joint].function[axis][dir].eq_count; 
for(  i  =  0;  i  <  eq_count;  i++ ) 

{ 

ier  =  fgets(  line,  132,  fp ); 

if(  ier  ==  NULL )  break;  /*  break  out  of  for  loop  7 
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sscanf(  line,  "%f  %e  %e  %e  %e", 

&table[joint]  .function[axis]  [dir]  .velocity  [i] , 
&tableOoint].function  [axis]  [dir]  .coeff[i]  [0] , 
&tableDoint].function[axis][dir].coeff[i][1], 
&tableDointj.function[axisj[dirj.coeff[i][2]1 
&tab le [jol n t] . fu n cti on [axi s] [dir], coeff  p] [3] ); 


} 

/*  break  out  of  the  while  loop  */ 
if( ier  ==  NULL)  break; 
axls++; 

table[joint].axis_count  =  axis; 

} 

fclose  (fp); 

}  /*end  for  each  joint*/ 


} 

int  InitJointCoeffTable  ( table ) 

/*iU***ftft  ***************  ***»*ft*0  *******  ft  *»*»«*«  **»*M**********^******A**********-»*' 

Function:  Given  a  pointer  to  the  coefficient  table,  this  function 
will  put  all  the  appropriate  initial  values  into  the  table. 

Parameters: 

table  -  pointer  to  the  coefficient  table  to  be  initialized. 

^********* ******* *****************  **  **********  *  *****************  *************************** 

CoeffTable  tableD; 

{ 

Int  I,  j,  k,  I,  m  ; 

/* for  each  joint  int  the  table*/ 
for  (m  =  0;  m  <  MAXJOINT;  m++) 

{ 

table[m].joint_name[0]  =  0; 


table[m].axis_count  =  0; 
table[m].axis[0]  =  0; 

table[m].axis[lj  =1; 

table[m].axis[2]  =  2; 


/* for  each  axis  in  that  joint*/ 
for(  I  =  0;  I  <  AXIS;  I++ ) 

{ 

/* for  each  direction  in  that  axis*/ 
for(j  =  0;j<DIR;j++) 

{ 

table[m].function[i][j].eq_count  =  0; 
table[m].function[i][j].direction_name[0]  =  0; 

/*for  each  velocity  in  that  direction*/ 
for(k  =  0;k<  VELOCITY;  k++) 

{ 

table[m].function[i][j].velocity[k]  =  0.0; 
for(  I  =  0;  I  <  COEFF;  I++ ) 

table[m].function[i][j].coeff[k][l]  =  0.0; 

} 
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UstCoeffTable(  table,  label,  count ) 

CoeffTable  tableQ; 
char  labelQ; 
int  count; 

/ . ************** - - * . . . . . ~***~~*~*~~*** 

Function:  Given  a  pointer  to  a  table,  this  function 
writes  the  function  to  stdout. 

Parameters: 

table  -  pointer  to  a  table. 

{ 

int  axis; 
int  dir; 

int  eq_count; 
int  i,j; 

printf  ( "%s\n“,  label); 
fa  (j  =  0 ;  j  <  count;  j++) 

{ 

printf(  "Joint  Name:  %s\n",  tableD].joint_name ); 
printf(  "DOFs  :  %d\n",  tabled-axis_count  )'< 
for(  axis  =  0;  axis  <  table[j].axis_count;  axis++ ) 

{ 

fa(  dir  =  0;  dir  <  2;  dir++  ) 

{ 

printf( "  Direction:  %s\n", 

tabled  function  [axis]  [dir]  .direction_name ) ; 

printf( "  Number  of  Functions:  %d\n", 
tabled.function[axis][dir].eq_count ); 

printf("  Velocity - Coefficients - \n" ); 

eq_count  =  tabled-function[axis][dir].eq_count; 
for(  i  =  0;  i  <  eq_count;  i++ ) 

{ 

printf("  %f  %E%E%E%E\n", 
tabled  function  [axis] [dir]  .velocity  p] , 
tabled  function  [axis]  [dir]  .coeffd  [0] , 
tabled-function[axis][dir].coeff[i][1], 
tabled  function  [axis]  [dir]  .coeffp]  [2] , 
tabled  function  [axis]  [dir]  .coeffp]  [3] ); 

} 

} 

} 

} 

} 

float  ipower(  x,  y ) 
float  x; 
int  y; 

l+*  +  kk  A*******  A*  kkk*i,kkk*Akiikk  ************************  AHA*  A*  *******  A****************  *****************  At  A  A 

float  ipower(  x,  y )  -  compute  x  to  the  y  power 
Function: 

Compute  the  value  of  x  raised  to  the  ylh  power. 

/*»»************»****»»*»»» *»**«*»*******»**»»« ****************  A*********** A* ****** *  +  ******  +  *+++**+++++++ 

{ 

float  z; 

2=1-0; 

while(  y- )  { z  *=  x; } 
retum(  z ); 
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MultiplyCoeff  (  coeffl  ,  coeff2,  value) 
float  coeffl  □; 
float  coeff2[|; 
float  value; 

J** *******  ************  **************************************** ************ ******************** AAA  *********** ******** 

Function: 

Multiply  the  array  of  coeffs  by  the  value  specified, 
coeffl  =  coeff2  *  value; 

Parameters: 

coeff  -  an  array  of  float  values  that  are  the  coefficients  to 

an  n  degree  polynomial. 

value  -  the  value  to  multiply  with. 

***********************************  *■*  ***************  *******************************************  *•*•■*  ***  *******  ******* 

{ 

int  i; 

for(  i  =  0;  i  <  COEFF;  i++  ) 

{ 

coeffl  [i]  =coeff2[i]  *  value; 

} 

} 

float  ComputeCoeff  ( coeff,  value ) 
float  coeff 0; 
float  value; 

I******* *****************  *******  ******* * ****************************** **********************************  ************ 

Function: 

Solve  the  coefficients  for  the  value  specified. 

Parameters: 

coeff  -  an  array  of  float  values  that  are  the  coefficients  to 
an  n  degree  polynomial. 

value  -  the  x  value  of  the  function  to  compute  the  y. 

**********  *******-************************************»*♦»************★******* ****************** A******************* 

{ 

int  i; 

float  answer; 

answer  =  coeff  [0]; 

for(  i  =  1 ;  i  <  COEFF;  i++) 

{ 

answer  +=  ipower(  value ,  i )  *  coeff p]; 

} 

return  (  answer ); 

} 

int  FTGetAxisNumber(  name ) 
char  name; 

I** ******************** A***************-*****-******************************************  **************  ****** ********** 

int  FTGetAxisNumber(  name )  Function: 

Given  an  axis  name  return  axis  number  associated  with  name. 
************************************************************************************************************************** 

{ 

char  a; 

a  =  toupper(  name ); 
if  ( a  ==  'X1 )  retum(  0 ); 
if  ( a  ==  Y )  retum(  1 ); 
if(a==’Z)  retum(2); 
retum(  -1 ); 
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} 

int  GetAxisDirection(  name ) 


Function:  Given  the  name  of  the  direction,  return  a  direction  number. 

Parameters:  name  -  direction  name. 

******************************************************************************  dHHHHHHHHHHHHHHHHHHHHHt  *  *+++++**  ****** 

char  ‘name; 

{ 

inti; 
int  k; 
i  =  0; 

while(  AxisMap[i].input_name[0] ) 

{ 

k  =  strlen(  AxisMap[i].output_name ); 

if(  stmcmp(  name,  AxisMap[i].output_name,  k )  ==  0 ) 

{ 

if(  AxisMap[i].direction  <  0 )  return(  AxisMap[i], direction  +  1  ); 
else  retum(  AxisMap[i].direction ); 

} 

i++; 

} 

return(2); 

} 

. . 

GetVelocityCode(  joint,  axis,  dir,  velocity ) 

Function: 

Get  the  code  number  associated  with  a  given  velocity 
found  in  the  force  function  velocity  table.  The  code 
number  is  selected  by  comparing  the  given  velocity  with 
the  velocities  in  the  table  and  chosing  the  index  number 
of  the  closest  value. 

int  GetVelocityCode(  force_table,  joint,  axis,  dir,  velocity ) 

CoeffTable  force_table[|; 

int  joint; 

int  axis; 

int  dir; 

float  velocity; 

{ 

inti; 

int  count; 
float  mid; 

/* 

size  of  the  velocity  is  the  same  as  the  number  of 
equations  within  the  given  axis  and  direction 

*/ 

count  =  force_table[joint].function[axis][dir].eq_count; 
f*  check  for  special  case  of  zero  vel.  */ 
f* no  equations*/ 
if  (count  ==  0 )  return  ( -1  ); 

if(  velocity  <  force_table[joint].function[axis][dir].velocity[0] ) 
retum(O); 

/*  scan  table  7 
for(  i  =  0;  i  <  count  - 1 ;  i++ ) 

{ 

if(  (velocity  >  force_table[joint].function[axis][dir].velocity[i])  && 
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{ 


(velocity  <  force_table(joint].function[axis][dir].velocityp+1]) ) 


/*  mid  =  (v[i+1]-v[i])/2.0  +  v[i]  */ 

mid  =  (force_table[joint].function[axis][dir].velocityp+1]  - 
force_tableOoint].function[axis][dir].velodty[i])  /  2.0  + 
force_tableOolnt].function[axis][dir].velocity[i]; 
f*  velocity  is  on  one  side  or  the  other  of  mid  point  */ 
if(  velocity  <  mid )  retum(  i ); 
retum(i  + 1 ); 

} 

} 

return ( count  - 1 );  /*  need  some  default  */ 

} 

int  GetCoeff(  joint,  axis,  dir,  velocity,  coeff ) 

Function: 

return  the  force  coefficents  associated  with  the  given 
joint,  axis,  direction  and  velocity. 

jirk  *■*  *******************************  **************  irk  *  *  *  *■■*■**■****'*'*'***1*  ******** 

int  GetCoeff(force_table,  joint,  axis,  dir,  velocity,  coeff ) 

CoeffTable  force_table[j; 

int  joint; 

int  axis; 

int  dir; 

float  velocity; 

float  coeffQ; 

{ 

int  vid; 
int  i; 

int  GetVelocityCodeO; 

vid  =  GetVelocityCode(  joint,  axis,  dir,  velocity ); 
for(  I  =  0;  i  <  COEFF;  i++) 

{ 

coeff  [i]  =  force_table[jolnt].function[axls][dir].coeff[vld][i]; 

} 

return(vid); 

} 

^****  ********  *******»★★*■*■* *************************  ******************************************************** **^ 
#include  <stdio.h> 
include  <math.h> 

#include  "lean.h" 
menu  0 

^****** *************  A********************************************************************************** ************* 

Function:  Print  a  menu. 

Parameters:  none. 

A*******************************************************************  ************************** kk  »**»»»»»  *»*»»★*★★*» 

{ 

printf  ( "  \n\n"); 

printf  ( "Fat  Menu  ===========>  :-)  \n"); 

printf  ( "0  ->  exit.  \n’); 

printf  ( "1  ->  Input  a  new  bf  and  weighAn\n“); 

printf  ( “2  ->  Print  Population  Coefficients^"); 

printf  ( "3  ->  Print  Lean  Body  Mass  coefficients^!"); 

printf  ( "4  ->  Print  This  individuals  coefficients\n\n"); 

printf  ( "5  ->  Write  out  (*.ffc  files)  individuals  coefficients^"); 

printf  ( "6  ->  Print  out  a  torque  summary  for  this  individual.\n"); 
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} 

main  (argc,  argv) 
int  argc; 
char*argvO; 

Function:  main  program.  Menu  selection  and  execution. 

Parameters:  body  fat  and  weight. 

***r*****rfr**fr*  *********************** ************************** ***************** 


Jaa ************* * ******** **** ***************** ********************** 

type  definitions  for  CoeffTable  data  structure  see  lean.h. 

Decription  of  force  table  data  structure. 

tabie(joint].function[axis][dir] 
force  functions  associated  with  a  joint,  these 
are  indexed  by  the  axis  number  and  the  direction 
code. 


tabie[joint]  .joint_name 
name  of  joint  80  characters  max. 
table[joint].axis[axis] 

tabie  of  axis  numbers  (DOFs)  for  the  joint. 
tabie[joint].axis_count 
number  of  axes  or  DOFs  for  the  joint. 
function[axis][dir].coeff[vid][coeff] 
force  function  coeffcients  index  by  the  velocity  code  (vid) 
(see  velocity  code)  and  coefficient  index  (1-4). 
function[axis][dir].velocity[vid] 
velocity  vaiue  associated  with  a  given  velocity  code  (vid). 
function[axis][dir].eq_count 
number  of  force  function  equations  associated  with 
a  given  axis  and  direction. 
function[axis][dir].direction_name 
name  of  a  given  direction,  i.e.  extension,  flexion,  etc. 


— - * . - . ********** - * - ** - 

CoeffTabie  pop_tabie[MAXJOINTj;  /*normaiized  popuiation  coefficients*/ 
CoeffTabie  lean_table[MAXJOINTj;  ^relationship:  mean  torque  to  lean  mass*/ 
CoeffTable  ind_tabie[MAXJOINT|;  /*Calculated  coeff  tabie  based  on 
weight  and  %body  fat*/ 

float  bf;  I*  body  fat  percentage*/ 

float  weight;  /*weight*/ 

int  done;  /loop  control  7 

int  option;  /*user  input  */ 

#if  (PC  ||  UNIX) 

/•process  command  line  arguments*/ 
if  ( argc  <  3 ) 

{ 


printf  ( "usage:  iean  <%body_fat>  <weight>  \n“); 
exit(  0 ); 


} 

#endif 
#if  MAC 


printf  ( "\n\n“); 

printf  ("input  body  fat  (0-1)  and  <weight(kg):"); 
scanf  ( “%f  %f,  &bf,  &weight); 
printf  ( “%f  %f ,  bf,  weight); 

#else 

bf  =  atof  ( argvjl] ); 
weight  =(float)  atof  ( argv[2j); 


"*★*■* 


7 


/ 
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#endif 
if  (bf>  1) 

{ 

printf  ( "  invalid  %  body  fat :-)  enter  value  between  0-1  \n"); 
exit(  0 ); 

} 

/*read  the  two  tables  1 .  Population  table  2.  lean  body  table*/ 
ReadCoeffTable  ( pop_table  ,  "ntc" ,  JC  ); 

ReadCoeffTable  ( leanjable,  "Ibc",  JC  ); 

/*calculate  this  individual’s  table*/ 

IniUointCoeffTable  ( ind_table ); 

CalculatelndividualCoeffTable 
(ind_table,  pop_table,  leanjable, 
bf,  weight,  JC ); 
done  =  FALSE; 
while  ( Idone ) 

{ 

menu  0; 

scanf  ( "%d",  &option); 
switch  (option) 

{ 

case  0:  /*exit */ 

printf  ("May  the  V’lean  force\"  be  with  you.  \n"); 
done  =  TRUE; 
break; 

case  1 :  /*lnput  new  values*/ 

printf  ("Input  body  fat  (0-1)  and  weight:"); 
scanf  ("%f  %f,  &bf,  &  weight); 

CalculatelndividualCoeffTable 
(ind_table,  pop_table,  leanjable, 
bf,  weight  JC); 
break; 

case  2:  /"Print  Population  Coefficients*/ 

UstCoeffTable  ( popjable,  "Population  Table",  JC); 

break; 

case  3:  /"Print  Lean  Body  mass  coefficients  7 

UstCoeffTable  ( leanjable,  "Lean  Body  mass  ->  average  torque",  JC); 
break;; 

case  4:  /"Print  Individuals  coefficients*/ 

UstCoeffTable  ( indjable,  "Individuals  coefficients",  JC); 
break; 

case  5:  /"Write  out  individuals  coefficient  file*/ 

WriteCoeffT ableT oFile(  indjable,  “ffc",  JC); 
printf  ( "  *.ffc  files  are  wntten\n"); 
break; 

case  6:  /"Write  out  a  summary  table  for  this  individual*/ 

TorqueSummary  ( leanjable,  bf,  weight,  JC); 
break; 
default; 

printf  ("  TRY  AGAIN,  (enter  options  0  ->6)  \n“); 
break; 

} 

} 

} 
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